Ground-state properties of one-dimensional anyon gases 
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We investigate the ground state of the one-dimensional interacting anyonic system based on the 
exact Bethe ansatz solution for arbitrary coupling constant (0 < c < oo) and statistics parameter 
(0 < K < tt). It is shown that the density of state in quasi- momentum k space and the ground 
state energy are determined by the renormalized coupling constant c'. The effect induced by the 
statistics parameter k exhibits in the momentum distribution in two aspects: Besides the effect of 
renormalized coupling, the anyonic statistics results in the nonsymmetric momentum distribution 
when the statistics parameter k deviates from (Bose statistics) and n (Fermi statistics) for any 
coupling constant c. The momentum distribution evolves from a Bose distribution to a Fermi one 
Eis K varies from to vr. The asymmetric momentum distribution comes from the contribution of 
the imaginary part of the non-diagonal element of reduced density matrix, which is an odd function 
of K. The peak at positive momentum will shift to negative momentum if k is negative. 
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I. INTRODUCTION 

The physical systems with fractional statistics have 
been a sul3ject of great interest [J, H, Although ini- 
tial studies of fractional statistics have focused on two- 
dimensional systems |1, 0, H, @, anyons (the particles 
obeying the fractional statistics) have also found applica- 
tion in various one-dimensional (ID) systems 0, 0,11, 0- 
Recently, many theoretical works have been dedicated to 
study the ID anyon gas 0, HI EH, El, [H [H, [H, [H, [13, 
[Tsl . [l§] since the integrable model of anyon gas_with a 5- 
function interaction was introduced by Kundu &!] despite 
no proof of the realization of anyon gas in ID. Because 
the statistical properties is related with the topological 
order, currently anyons also draw intensive attention in 
topological quantum computation. Particularly, with the 
rapid progress in the regime of cold atom and with the 
excellent controllability of the neutral atoms trapped in 
the optical lattice, cold atom has become a popular ex- 
perimental platform in many research regimes. Rotating 
Bose-Einstein condensates and cold atoms in optical lat- 
tice have also been proposed to create, manipulate, and 
test anyons jH, HH, ill] . 

On the other hand, the experimental progress of 
trapped ID cold atom systems [2^ [H, [H, HBj has trig- 
gered more and more theorists to study the many body 
physics of ID correlate d sy stems beyond the mean-field 
theory [13, [H, H [M (H. Further, the ability of 
tuning the effective ID interactions by Feshbach reso- 
nance leads experiment accessible to the strong correla- 
tion regime [25L [26l]. In the light of recent experiments, 
studies of the integrable model have re-attracted much 
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attention today [H, [s^, [H, [H, [13] because with its exact 
solution we can obtain the properties in the full physical 
regime even in the strong correlated limit. By apply- 
ing the exact solution of Lieb-Liniger to the ID inter- 
acting Bose gas in a hard- wall trap [H, [s^ , the evolu- 
tion from the Bose-Einstein condensate to "fermionized" 
Tonks-Girardeau (TG) gas was explicitly and exactly dis- 
played in our previous works . 35. , .36] , which has also been 
confirmed by various numerical methods [13, [4l[, . Al- 
though the density distribution of Bose gas in the limit 
of infinite repulsion is identical to the free Fermi gas [i^ , 
the momentum distributions exhibit quite different prop- 
erties from the free Fermi gas because of their respective 
statistics [3 [4^ . It is well known that the different per- 
mutation symmetry leads to the physics of Bose gases in 
many aspects differing from the physics of Fermi gases. 
As a natural generalization, the anyon gas interpolates 
between the Bose and Fermi gas. No doubt the inte- 
grable anyon model provides us the possibility to study 
the effect of generalized permutation symmetry in an ex- 
act manner although it is still a theoretically hypothetical 
model. So far much attention has focused on its ther- 
modynamic properties and correlation functions [lol . [Tl| , 
however no general result for the momentum distribution 
has been given except for the anyon gas in the infinite in- 
teraction limit (or the anyonic TG gas) [H, [11] . In this 
paper, we investigate the ground state of the ID any- 
onic system with finite number of anyons in the whole 
regime of interaction parameter, and particularly we shall 
study how the fractional statistics affect the ground state 
properties, such as the ground state energy and the mo- 
mentum distribution. By numerically solving the Bethe 
ansatz equations, we obtain the ground state wave func- 
tion and thus the momentum distribution for different 
coupling constant and statistics parameter. The obvi- 
ous properties of anyonic statistics are displayed in the 
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momentum distribution and with the change of statistics 
parameter the system exhibits the Bose statistics, Fermi 
statistics and the fractional statistics in between. 

The paper is organized as follows. In Sec. II, we give 
a brief review of ID anyonic model and introduce its so- 
lution with Bethe ansatz method. In Sec. Ill, we present 
the density of ground state in quasi-momentum space, 
ground state energy and momentum distributions for dif- 
ferent coupling constant and statistics parameter. A brief 
summary is given in Sec. IV. 



II. FORMULATION OF THE MODEL AND ITS 
SOLUTION 



We consider the one dimensional system of N anyons 
in length L with second quantized Hamiltonian 



Ha = 



2m 



dx 



dx dx 



+ ^ £ dx^\ (X) (X) ^A (X) ^A (X) , (1) 

in which m is the mass of anyon and giD denotes the 
effective interacting constant between anyons. The field 
operator Vl/^ {x) and '^a (x) obey the anyonic commuta- 
tion relations 

^Aixi)^\{x2) = e-''''^^^-^'^^\{x2)^A{xi) 

+S (xi - X2) , 

where the sign function e (x) gives —1, or 1 depending 
on whether x is negative, zero, or positive, k is the pa- 
rameter related with statistics and will be restricted in 
the regime [0,7r] in the present work. Particularly, k = 
and TT correspond to Bose statistics and Fermi statistics, 
respectively. A standard rescaling procedure brings the 
Hamiltonian into a dimensionless one 



Ha 



dx 



dx dx 



c*^ {x) (x) -^A {x) -^A [x) 



with c being the dimensionless coupling constant. The 
Hamiltonian take the same form as the Lieb-Liniger 
Bose model except that the field operator satisfy anyonic 
statistics. Similar to the Bose model the eigen problem 
for the Hamiltonian can be reduced to solve the problem 
of A^-anyons with (5-interaction [1, 



H^{xi,--- ,xn) =£;*(xi,--- ,xn) 



with 



N 



92 



i=i J 3<i 



Xl] 



(2) 



(3) 



The anyonic commutation relations require the wavefunc- 
tion satisfy the generalized symmetry 



* (xi, • • • ,x-,-, • • • ,X/, • • • ,Xn) 
= e""^*(xi,--- ,x/,--- ,Xj,--- ,xn) 

with the anyonic phase 



(4) 



H ({xj-xk)- ^ 

k=j+l k=j+l 



e{xi - Xk) 



Using the coordinate Bethe ansatz we can obtain the 
wavefunction 



4'(xi, 



,xn) = ^( 



X exp (x 

X (pQ (Xgj , Xq2 , • • • , Xqjy ) . 



qi 1 -^92 J 



' (2^132 Xq-^ ) 

(5) 



where Q is used to label the region < < Xq^ < 
■ ■ ■ < Xq^ < L and 6{x — y) is the step function. In 
different region the phase factors of the wavefunction are 
determined by A = X]i<j ^ (^i ~ which results in the 
exchange symmetry dictated by eq. The wave func- 
tion ifQ (xi, ■ • ■ , xjv) is taken as the Bethe ansatz type 



IPQ [Xql, ■ ■ ■ ,XqN) = ^ 



Ap^...p^ exp i^fc; 



(6) 

In order to get a physical result, the twisted boundary 
condition 



^'(0,X2, 



,Xn 



)^e'-(^~''^^{L,X2,---XN) 



is used here. Under the twisted boundary condition, we 
have the Bethe ansatz equations 



exp (ikjL) 



N 

n 



iki — ikj + c' 



(7) 



with j ~ 1,2, ■ ■ ■ , N, where 

c' — c/ cos (k/2) 

is the renormalized coupling constant and k the statistics 
parameter. It is deserved to notice that the system is 
reduced to the Lieb-Liniger Boson model when k = 
and the system is reduced to the free Fermi one when 
K — TT. The coefficients Ap take the form of 



N 

=ep Yl {ik, 



Pi ikp. -\- c 



in which Gp denotes a -!-(—) sign factor associated with 
even (odd) permutations. With the set of quantum num- 
ber {fcj } known as quasimomentum the energy eigenvalue 
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and the total momentum can be formulated as 

N 



(8) 



and P — Taking the logarithm of Bethe ansatz 

equations, we have 



N 



kjL — HjTT— farctan^^^ — j ' 



(9) 



In the following only the case of c > will be investigated. 
The set of integer {rij} determine an eigenstate and for 
the ground state Uj ^ j ~ {N + 1) /2 {1 < j < N). 
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III. PROPERTIES OF THE GROUND STATE 

By numerically solving the set of transcendental equa- 
tions eq.® , the quasimomentum {kj} and thus the wave 
function can be decided exactly. For c > 0, the set of ki is 
unique and real, therefore the density of state in k space 
p{k) can be formulated as [1^ 

In Fig. 1, we plot the density of ground state in quasi- 
momentum k space with the interaction parameter c fixed 
to c = 0.1 for several k. It is shown that the density 
of state of anyons with coupling constant c and statis- 
tical parameter k is identical to the density of state of 
Lieb-Liniger Bose model with effective coupling constant 
c' = c/cos{k/2). This is clear by comparing the Bethe 
ansatz equation ^ with that of Lieb-Liniger model [s^ . 
When K = 0, the model is reduced to Bose model and k 
mainly distribute in the regime around zero. While for 
the case of k = tt, we have kj = (j — {N + 1)/2)2tt/L 
with j ~ 1, • • • , N even for the weak coupling c — 0.1, 
which is the same as the situation of c' = oo of the Lieb- 
Liniger Bose model. Here the statistics parameter k af- 
fects the density of state in quasi-momentum space only 
by renormalizing the coupling constant c into c'. For 
the case of c = 0.1, we note that when k approaches tt 
minor change of k will result in the great change of the 
density of ground state in k space. For instance, there 
are obvious difference among the densities corresponding 
to K — O.Ott, 0.997r, and tt due to the dramatic changes 
of the effective coupling constant c' = 0.64, 6.37 and oo 
correspondingly. In Fig. 2, we display the ground-state 
energy versus k for different coupling constant c. The 
ground-state energy increases with the increase of k and 
the coupling constant c. For the weak coupling the energy 
of the system is small for most k but it abruptly increases 
and arrives at the maximum when k approach tt. In the 
strong coupling limit (c — s- oo) the ground-state energy 
is almost a constant independent of the statistics param- 
eter. Similarly the ground-state energy of the anyon gas 



FIG. 1: (color online) The density of state in quasi- 
momentum k space for the ground state for N = 200, L — 200 
and c = 0.1. k = (solid line), k = O.Stt (dashed hnes), 
K — O.Ott (dot hnes), k — 0.997r (dash dot hnes) and n — n 
(dash dot dot lines). 
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FIG. 2: (color online) The energy of ground state for A'' — 200 
and L — 200. c = 0.1 (solid line), c = 1 (dashed lines), c — 5 
(dot hnes), c = 10 (dash dot hnes), c = 50 (dash dot dot 
hnes), c — 1000 (short dash hnes) and c — oo (short dot 
hnes) . 



with coupling constant c and statistical parameter k is 
identical to that of the corresponding Bose gas with ef- 
fective coupling constant c' = c/cos{k/2). 

Now it is clear that the density of state in k space 
and the ground-state energy are determined completely 
by the renormalized coupling constant c'. Furthermore, 
from the eqs. ([U and ([9]), we see that the energy level 
structure of the anyon gas is exactly the same as the 
corresponding bosonic model with renormalized coupling 
constant c'. Therefore the thermodynamic properties of 
the anyon gas are the same as the well known thermody- 
namic properties of the ID boson gas [H, with the 
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FIG. 3: (color online) Momentum distribution of anyons for 
iV = 4, L = 1 and c = 0.1. k = (solid line), k = Stt/S 
(dashed lines), k = Stt/S (dot lines), k — 7tt/8 (dash dot 
lines), K — 3l7r/32 (short dash lines) and k — n (short dot 
lines) . 



FIG. 4: (color online) Momentum distribution of Bosons for 
iV = 4, L = 1 and c = 0.1/cos(k/2). k = (solid line), 
K = 37r/8 (dashed lines), k — 5n/8 (dot lines), k = Ttt/S 
(dash dot lines), k = 3l7r/32 (short dash lines) and k = n 
(short dot lines). 



effective coupling c' = c/cos{k/2). This implies that the 
intrinsic difference between the anyon gas and the cor- 
responding Bose gas is hard to be distinguished just by 
the thermodynamic properties. However, due to the dif- 
ferent exchange symmetry of the wave functions, the ob- 
servables associated with the wave functions rather than 
the square of wave functions, such as the off-diagonal re- 
duced density matrix and the momentum distributions, 
should display quite different behaviors. In terms of the 
ground state wave function \1/ (xi, • • • , xn), the one-body 
reduced density matrix is given by 

p{x,x') 

_ N dx2 ■ ■ ■ dxN^* {x,X2, • • • , a;Ar) 5' {x' ,X2, • • • , xn) 

dxi ■ ■ -dxN I* (xi,a;2, • ■ • ,XN)f' 
while the momentum distribution 

n{p)^—[ dx [ dx'p{x,x')e~'P('=-''') (11) 
27r Jo Jo 

is the Fourier transformation of p{x, x'). Below we rescale 
the momentum distribution n{p) and momentum p into 
dimensionless forms in the units of L/27r and 27r/_L re- 
spectively and the original notations are reserved. 

To give a concrete example, we display the momen- 
tum distribution of the anyon gas with coupling constant 
c = 0.1 and TV = 4 in Fig. 3 for various statistics parame- 
ters. At K = 0, the momentum distribution shows a typi- 
cal Boson-type distribution with an obvious peak around 
the zero momentum point, whereas the momentum dis- 
tribution is identical to that of free spinless Fermions at 
K = TT, exhibiting shell structure. In these two situa- 
tions the momentum distributions are symmetric. When 
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FIG. 5: (color online) Momentum distribution of anyons for 
= 4, L = 1 and c — oo. n — (solid line), k = tt/4 (dashed 
lines), K, — Tr/2 (dot lines), k — Sn/A (dash dot lines) and 
K = TT (short dash lines). 



K deviates from these two extrema points, the momen- 
tum distribution becomes asymmetric: the anyons dis- 
tribute with more probability in the regime of positive 
momentum whereas with relatively small probability in 
the other side. With the increase of k, the distribution 
redistributes between the regime of positive and negative 
p, and gets wider with more and more peaks appearing. 
Finally, N peaks appear at k = tt and the momentum 
distribution shows the characteristic shell structure of 
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Fermions. 

As a comparison, in Fig. 4 we display the momentum 
distribution of the Lieb-Liniger boson gas for = 4 with 
the renormahzed coupHng constants c = 0.1/cos(k/2) = 
0.1,0.12,0.18,0.51,2.04 and oo corresponding to those 
for anyons in Fig. 3, whore we have c = 0.1 and 
K = 0, 37r/8, 57r/8, 77r/8, 3l7r/32, and tt respectively. The 
momentum distributions for c = 0.1, 0.12, 0.18, 0.51, 2.04 
are not obviously different in Fig. 4 as they all fall in the 
weakly interacting regime. We also display the momen- 
tum distribution for c = 50 in Fig. 4. The height of the 
peak shrinks with the increase of the repulsive interaction 
and stronger repulsive interaction between the bosonic 
atoms tends to spread out the distribution widely. Even 
in the TG limit, the momentum distribution of bosons 
docs not show shell structure like free fermions. For all 
coupling constant, bosons always congregate symmetri- 
cally in the regime around p = with the most probabil- 
ity and with the increase of momentum p the probability 
diminishes rapidly. Obviously, the asymmetry in the mo- 
mentum distribution of anyons can not be attributed to 
renormalization of the effective coupling constant and is 
induced completely by the the permutation symmetry re- 
lated to statistic parameter k, which can be regarded as 
the indication of anyons. 

In order to clarify more clearly the effect induced by 
the statistic parameter k, we show the momentum dis- 
tributions for the system with N = 4 and c = oo in 
Fig. 5. In this case the renormalized coupling constant 
c' = 00 and the quasimomentum {kj} shall be same for 
arbitrary k, therefore the statistics parameter only re- 
sults in different permutation symmetry embodied in the 
wavefunction. There appears the similar situation as the 
case for c = 0.1 in Fig. 3. In the Bosc limit (k = 0) and 
in the Fermi limit (k = tt) the momentum distributions 
are symmetric while in between the distributions take on 
the nonsymmetric profiles and evolve from a Bose distri- 
bution into a free spinless Fermi distribution with the in- 
crease of K. The difference from the case of weak coupling 
is that the system distributes in a wider regime relative 
to that of weak coupling. Mathematically, the asymmet- 
ric momentum distribution is induced by the existence of 
imaginary part of the density matrix (see appendix for 
detail), which is an odd function of statistical parameter 
such that the peak at positive momentum in Fig. 3 and 
Fig. 5 is a result of the choice of positive k, and the peak 
will shift to negative momentum if k is negative. 



IV. CONCLUSIONS 

In summary, we have investigated the ground-state 
properties of ID anyon gas based on the exact Bethe 
ansatz solution for arbitrary coupling constant (0 < c < 
oo) and statistics parameter (0 < k < tt). The density 
of state in quasimomentum space and the ground-state 
energy arc determined by the renormalized coupling con- 
stant c'. The anyonic system with coupling constant c 



and statistics parameter k has the same density of state 
and energy as those of the Lieb-Liniger Bose model with 
coupling c' = c/cos{k/2). Besides the renormalization 
of the effective coupling constant, the additional effect 
induced by the statistics exhibits in the momentum dis- 
tribution. While in the limit of k = and k = it, the 
anyon gas is reduced to the Bose gas and free spinless 
Fermi gas, respectively, with both of momentum distri- 
butions being symmetric. When the statistics parameter 
deviates from these two points, the momentum distribu- 
tion is asymmetric and evolves from the Bose distribution 
to Fermi distribution with the increase of k. Although 
the thermodynamic properties of anyon sytem are even 
functions of the statistic parameter k because the renor- 
malized coupling constant is even function of it, the mo- 
mentum distribution is not so because of the existence of 
the imaginary part of the non-diagonal part of the one 
body density matrix, which is an odd function of k. If 
K is negative the momentum distribution is equal to the 
distribution for positive k with the mapping between the 
positive momentum and negative momentum. Our re- 
sults interpolate between the known results of Bose and 
Fermi gas and clarify the effect of generalized permuta- 
tion symmetry of the anyon gas in an exact manner. 



Acknowledgments 

S.C. would like to appreciate X. W. Guan for useful 
discussions. This work is supported by NSF of China 
under Grants No. 10574150 and No. 10774095, National 
Program for Basic Research of MOST China, the 973 
Program under Grant No. 2006CB921102, and Shanxi 
Province Youth Science Foundation under Grant No. 
20051001. 



APPENDIX A 

The asymmetry of the momentum distribution can be 

attributed to the special permutation symmetry of the 
anyon wavefunction characterized by the statistical pa- 
rameter K. To see it clearly, we display the case of two 
anyons in detail in this appendix. For the two anyon 
systems, the Bethe ansatz wavefunction takes the form 
of 

^ixi,X2) = 9{xi < X2)exp{iK/2)ip{xi,X2) + 
9 {x2 < xi) exp {—in/2) {x2, x\) , 

where ip{xi,X2) = A12 cxp(v'fcia;i -I- ik2X2) + 
y4.2i exp(zA:2-Ti + ikiX2). The reduced one body density 
matrix p{x,x') = (1x2^^* {x,X2)'^ {x',X2) can be 
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evaluated with wavefunction, i.e., 
p{x,x') 

= / dx2[0{x < X2)0{x' < X2)(p* {x,X2)(p{x',X2) 

Jo 

+6 {x < X2) 6 {x2 < x') exp {-in) tp* {x, X2) ip ix2,x') 
+9 {x2 < x) 9 {x' < X2) exp (in) (p* {x2,x) p {x' , X2) 

+9 {X2 < x)9 {X2 < x') p* {X2,x) (p {X2,x')] . 

When the statistical parameter k deviates from and 
TT, the density matrix is always complex for x 7^ 
x' . From the definition of reduced density matrix, we 
have p{x,x') — p*{x',x), which implies He[p{x, x')] = 
Re [p{x' , x)] and Im [p{x, x')] = — Im [p{x' , x)]. Therefore 



the momentum distribution can be represented as 

n{p) = ^ t dx dx'p{x,x')e-'P(''-^') 
27r Jo Jq 

= — [ dx [ dx' {Re[p{x,x')]cosp{x — x') 
27r Jo Jo 

+ Im [p{x, x')] sinp (x ^ x')} . 

The second term is the odd function of momentum p 
such that the momentum distribution becomes asymmet- 
ric about p. Furthermore, the imaginary part Im [p{x, x')] 
is an odd function of k. Therefore, if we take k as nega- 
tive (k —>■ — k), the peak at positive momentum as shown 
in Fig. 3 and Fig. 5 will shift to negative momentum. 
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